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The standard cosmological model, based on general relativity with an inflationary era, is very 
effective in accounting for a broad range of observed features of the universe. However, the ongoing 
puzzles about the nature of dark matter and dark energy, together with the problem of a fundamental 
theoretical framework for inflation, indicate that cosmology may be probing the limits of validity 
of general relativity. The early universe provides a testing ground for theories of gravity, since 
gravitational dynamics can lead to characteristic imprints on the CMB and other cosmological 
observations. Precision cosmology is in principle a means to constrain and possibly falsify candidate 
' quantum gravity theories like M theory. Generalized Randall-Sundrum brane-worlds provide a 

■ phenomenological means to test aspects of M theory. I outline the 1+3-covariant approach to 
' cosmological perturbations in these brane-worlds, and its application to CMB anisotropies. 

(N : 

Ph' I. INTRODUCTION 

<\ 

, M theory, the 1 l-dimensional theory that encompasses the known superstring theories, is only partially understood, 

■ but is widely considered to be a promising potential route to quantum gravity [1], and is therefore an important 
\ candidate for cosmological testing. Currently there are not realistic M theory cosmological solutions, so that it is 
• reasonable to use simplified phenomenological models that share some of the key features of M theory, especially 

^ ' branes. In brane cosmology, the observable universe is a l-|-3-dimensional "brane" surface moving in a higher- 
dimensional "bulk" spacetime. Fields and particles in the non-gravitational sector are confined to the brane, while 
gravity propagates in the bulk. The simplest, and yet sufficiently general, phenomenological brane-world models are 
the cosmological generalizations [2-4] of the Randall-Sundrum II model [5]. In the RSII brane-world, the bulk is 
5-dimensional anti-de Sitter spacetime, so that the extra dimension is infinite. The generalized brane-worlds (see [6,7] 
CO , for recent reviews) also have non-compact extra dimension.^ The other 6 extra spatial dimensions of M theory may 

■ be assumed to be stabilized and compactified on a very small scale, so that they do not affect the dynamics over the 
Q I range of validity of the brane-world model, i.e. for energies sufficiently below the string scale, for which the brane 

may be treated as infinitely thin. The RSII models have the additional advantage that they provide a framework for 
^ investigating aspects of holography and the AdS/CFT correspondence. 

bi). What prevents gravity from 'leaking' into the infinite extra dimension at low energies is the negative bulk cosmo- 

■ logical constant. 
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where £ is the curvature radius if the bulk is AdSs. Corrections to Newton's law in the weak-field static limit are 
0{e/r^) [5]: 



GM ( 2£' 



Experiments currently impose an upper bound £ ^ 1 mm. On the brane, the negative A5 is offset by the positive 
brane tension A. The effective cosmological constant on the brane is 



^Various extensions of the generalized RSII brane-worlds are not discussed here. The simplest extension is to introduce a 
second brane, so that the extra dimension is compact (but much larger than Planck scale), or a scalar field (gravitational 
sector) in the bulk, or both. See [7] for further discussion and references (including the "ekpyrotic" and cyclic models where the 
potential of the bulk scalar causes branes to collide, which may initiate a big bang and provide an alternative to inflation [8]). 
Other extensions involve corrections to the action, including 4-dimensional "induced gravity" and 5-dimensional Gauss-Bonnet 
corrections; see, e.g. [9]. 
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A=1(A5 + /^2A), ^ = 8nG=^, (3) 

where M4 ~ 10^^ GeV is the effective Planck scale on the brane. This is not the true fundamental gravity scale, which 
can be much lower, offering the possibility of a resolution of the hierarchy problem, as well as the exciting prospect 

that quantum gravity effects could be observable in particle accelerators and cosmic ray showers. The fundamental 
energy scale can be as low as ~ TeV in some brane-world scenarios, but in generalized RSII models it is higher, 
M5 > 10^ TeV, and is related to M4 via 

M| = ^. (4) 

The bound £ < 1 mm implies that A is above the electroweak scale, A > (100 GeV)^. At high energies {p » A) 
in the early universe, gravity becomes 5-dimensional and there are significant corrections to standard cosmological 
dynamics. There are also corrections that can operate at low energies, mediated by bulk graviton or Kaluza-Klein (KK) 
modes. Both types of correction play an important role in cosmological perturbations. In particular, 5-dimensional 
gravitational-wave modes introduce nonlocal effects from the viewpoint of brane-bound observers [10,11]. 



Brane-world inflation 



The unperturbed cosmological brane-world is a Friedmann brane in a Schwarzschild-AdSs bulk [2,12]. High-energy 
brane-world modifications to the dynamics of inflation on the brane have been investigated [4,13]. Essentially, the 
high-energy corrections provide increased Hubble damping. 



y(<p)>A => H 
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(5) 



thus making slow-roll inflation possible even for potentials V{(fi) that would be too steep in standard cosmology [4,14]. 
This can be seen clearly from the slow-roll parameters {V A) 
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where figi , r/gr arc the standard general relativity slow- roll parameters. Steep potentials can inflate at high energy and 
then naturally stop inflating when V drops below A. These models can be constrained because they typically generate 
a blue spectrum of gravitational waves which can disturb nucleosynthesis [14] . They also raise the intriguing prospect 
that the inflaton could act as dark matter or quintessence at low energies [14,15]. 

Large-scale scalar perturbations generated by slow-roll inflation {V ^ A) have an enhanced amplitude compared 
with the standard general relativity case [4]: 
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This means that COBE-scale perturbations can be generated when the inflaton is well below M4. For example, 
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V^cobe 
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M4 « M4 , 



for V = ^rn^ip^. The scalar spectral index is in general given by 

rig = 1 — 6e -|- 2ry , 



(7) 



(8) 



(9) 



and is driven closer to 1 (compared to general relativity) by high-energy effects. 

High-energy inflation on the brane also generates a zero-mode (4-dimensioiiaI graviton mode) of tensor perturba- 
tions, and stretches it to supcr-Hubble scales. This zero-mode has the same qualitative features as in general relativity, 
remaining frozen at constant amplitude while beyond the Hubble horizon. Its amplitude is enhanced at high energies, 
although the enhancement is much less than for scalar perturbations [16]: 



2 



_8V_ 

25M| 



(10) 
(11) 



Equation (11) means that brane- world effects suppress the large-scale tensor contribution to CMB anisotropics. The 
tensor spectral index has a smaller magnitude than in general relativity, but obeys the same consistency relation: 



(12) 



The massive KK modes (5-dimensional graviton modes which have an effective mass from a brane observer view- 
point) remain in the vacuum state during slow-roll inflation [16,17]. The evolution of the super-Hubble zero mode is 
the same as in general relativity, so that high-energy brane-world effects in the early universe serve only to rescale 
the amplitude. However, when the zero mode re-enters the Hubble horizon, massive KK modes can be excited. This 
may be a very small effect, but it remains to be properly quantified. 

Vector perturbations in the bulk metric can support vector metric perturbations on the brane, even in the absence of 
matter perturbations. However, there is no normalizable zero mode, and the massive KK modes stay in the vacuum 
state during brane-world inflation [18]. Therefore, as in general relativity, we can neglect vector perturbations in 
inflationary cosmology. 



Perturbation evolution 



The background dynamics of brane-world cosmology is known exactly. Large-scale cosmological perturbations on 
the brane are well understood [4,11,16,19,20]. However, without a solution for small-scale perturbations, we remain 
unable to predict the CMB anisotropics in brane-world cosmology, and the CMB provides the key means to test 
the scenario. The problem is that the 5-dimensional bulk perturbation equations must be solved in order to solve 
for perturbations on the brane. The 5-dimensional equations are partial differential equations for the 3-dimensional 
Fourier modes, with complicated boundary conditions. In fact, even the Sachs- Wolfe effect requires information from 
the 5-dimensional solutions; although the large-scale density perturbations can be determined without knowing the 
5-dimensional solutions [11,19], the Sachs-Wolfe effect requires the large-scale metric perturbations, and these are 
related to the density perturbations in a way that involves the KK modes [20] . 

The theory of gauge-invariant perturbations in brane-world cosmology has been extensively investigated and de- 
veloped [4,10,11,14,16-24] and is qualitatively well understood. The key remaining task is integration of the coupled 
brane-bulk perturbation equations; up to now, only special cases have been solved, where these equations effectively 
decouple. In general, and for the crucial case of calculating CMB anisotropics [20,22 24], the coupled system must be 
solved. From the brane viewpoint, the bulk effects, i.e. the high-energy corrections and the KK modes, act as source 
terms for the brane perturbation equations. At the same time, perturbations of matter on the brane can generate 
KK modes (i.e., emit 5-dimensional gravitons into the bulk) which propagate in the bulk and can interact with the 
brane. This nonlocal interaction amongst the perturbations is at the core of the complexity of the problem. It can 
be elegantly expressed via integro- differential equations [10], which take the form 

Ak{t) = y" dt'g{t,t')Bk{t'), (13) 

where G is the bulk retarded Green's function evaluated on the brane, and A). , are made up of brane metric and 
matter perturbations and their (brane) derivatives, and include high-energy corrections to the background dynamics. 



II. COVARIANT DYNAMICS AND PERTURBATIONS 



The 5D field equations are 



= -As (^Ws + 5{y) ^ [-XgAB + Tab] , (14) 



where y is a Gaussian normal coordinate orthogonal to the brane, which is at y = 0, the induced metric on {y = const} 
is Qab = ^^^dAB — uaUb with the unit normal, and Tab is the energy- momentum tensor of particles and fields 
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confined to the brane (with TAsn^ = 0). The effective field equations on the brane are derived from the Gauss-Codazzi 
equations and the Darmois-Israel junction conditions (using Z2-symmetry) [3]: 

Gab = -^9ab + K^Tab + 6y 5^6 - Sab , (15) 

where Sab ~ (Tab)'^ is the high-energy correction term, which is neghgiblc for p <C A, while 8ab is the projection 
of the bulk Weyl tensor on the brane. This term encodes corrections from KK or 5D graviton effects. From the 
brane-observer viewpoint, the energy-momentum corrections in Sab are local, whereas the KK corrections in Sat arc 
nonlocal, since they incorporate 5D gravity wave modes, as discussed above. These nonlocal corrections cannot be 
determined purely from data on the brane, and so the effective field equations are not a closed system. One needs to 
supplement them by 5D equations governing Sab, which are obtained from the 5D Einstein and Bianchi equations [3]. 

The trace free Sab contributes an effective energy density p* , pressure /9*/3, momentum density g* and anisotropic 
stress w*^ on the brane. In a 1+3-covariant decomposition, 



~^Sab = P* iuaUb + ^hab ] + Q^Ub + q^Ua + , (16) 



where u°' is a physically determined 4-velocity on the brane and hab = gab + UaUt projects into the comoving rest 
space at each event. The KK anisotropic stress n*^ incorporates the spin-0 ("Coulomb"), spin-1 (gravimagnetic) and 
spin-2 (gravitational wave) 4D modes of the 5D graviton. The KK momentum density q* incorporates spin-0 and 
spin-1 modes, and the KK energy density p* (the "dark radiation") incorporates the spin-0 mode. The brane "feels" 
the bulk gravitational field through these terms. In the background, q^ = = tt*^, since only the "dark radiation" 
term is compatible with Friedmann symmetry. 

The brane-world corrections can conveniently be consolidated into an effective total energy density, pressure, mo- 
mentum density and anisotropic stress. Linearizing the general nonlinear expressions [11] we obtain 



2A p 

+ + f ' 

qf =qa{l + £)+q:, (19) 
<t =^ab{l-'-^)+Kb, (20) 

where p = '^^p^^^ and p = X^jP'^*'' arc the total matter-radiation density and pressure, qa = ^iqa^ is the total 
matter-radiation momentum density, and i^ab is the photon anisotropic stress (neglecting that of neutrinos, baryons 
and CDM). 



Energy-momentum conservation, 



together with the 4D Bianchi identity, lead to 



VX6 = 0, (21) 



V'faft = — W-Sab , (22) 

which shows qualitatively how 1-1-3 spacetime variations in the matter-radiation on the brane can source KK modes. 
The 1+3-covariant decomposition of Eq. (21) leads to the standard energy and (linearized) momentum conservation 
equations, 

p + e(p + p)+D^„ = 0, (23) 

qa + ^Hqa + BaP + [p + p)Aa + B^'TTab = , (24) 

where Q is the volume expansion rate, which reduces to 3H in the background {H is the background Hubble rate), Aa 
is the 4-acceleration, and Do is the covariant derivative in the rest space (i.e. DaF*" "...c = ha'^h!'e ■ ■ ■ hJV dF^'" ■■■})■ 
The absence of bulk source terms in the conservation equations is a consequence of having A5 as the only 5D source 
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in the bulk. If there is a bulk scalar field, then there is energy-momentum exchange between the brane and bulk (in 
addition to the gravitational interaction) [25]. 

Equation (22) may be thought of as the "nonlocal conservation equation" . Linearizing the general 1+3-covariant 
decomposition [11], we obtain 

p* + lep*+DX=0, (25) 



3 



-D„p + 3i?ga + ^D''7r„b 



(26) 



At linear order, spatial inhomogeneity (D^p), peculiar motions {qa = pVa) and anisotropic stresses (iTab) in the matter- 
radiation on the brane are seen to be sources for KK modes (or 5D graviton emission into the bulk) . Qualitatively and 
geometrically this can be understood as follows: the non-uniform 5D gravitational field generated by inhomogcneous 
and anisotropic 4D matter-radiation contributes to the 5D Weyl tensor, which nonlocally "backreacts" on the brane 
via its projection Sab- Note also that the source terms are suppressed at low energies, and during quasi-de Sitter 
inflation on the brane. 

Equations (25) and (26) arc propagation equations for p* and q*. There is no propagation equation on the brane 
for TT*^; if there were such an equation, then one could determine the KK modes purely from data on the brane, which 
would violate causality for 5D gravitational waves. 

In the background, the modified Priedmann equations are 

H' = ^P^' +Ia+^, (27) 
3 3 a"^ 



H = --{p^«+p^^) + ^. (28) 
By Eq. (25), the KK energy density behaves like dark radiation: 

p* (X 1 . (29) 



The source of the dark radiation is the tidal (Coulomb) effect of a 5D black hole in the bulk. When the black hole mass 

vanishes, the bulk geometry reduces to AdS^ and p* = 0. In order to avoid a naked singularity, we assume that the 
black hole mass is non-negative, so that p* > 0. This additional effective relativistic degree of freedom is constrained 
by nucleosynthesis and CMB observations to be no more than -^3% of the radiation energy density [20,22,26]: 



Prad 



< 0.03 (30) 



If p* = and K = = A, then the exact solution of the Priedmann equations is [2] 

a = const [t{t + tA)]'/^^"'+'^ , tx = ^ < 10"^ sec , (31) 

vttA 



where w = p/p is assumed constant. If p* 7^ (but K = = A), then the solution for the radiation era (w = |) is [22] 

a=const[t{t + tx)]'^\ tx= ^ rr^^'w (32) 

AVttX (1 -I- p*/p) 



For t ^ t\ we recover from Eqs. (31) and (32) the standard behaviour, a oc ^2/3(^-1-1)^ whereas for t <C tx, we have 
the very different behaviour, a oc 

In the 1+3-covariant description of perturbations [11], wc isolate the KK anisotropic stress Tr*^ as the term that 
must be determined from 5D equations. Once tt*;, is determined in this way, the 1-1-3 perturbation equations on the 
brane form a closed system. The KK terms act as source terms modifying the standard general relativity perturbation 
equations, together with the high-energy corrections. For example, the propagation equation for the shear is [11] 

aab + 2Haab + Eab - yT^afe - D(oA) = -^Kb ~ + 3"^))^ '^"■b , (33) 
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where Eab is the electric part of the 4D brane Weyl tensor (not to be confused with £ab)- In general relativity, the 
right hand side is zero. In the brane-world, the first source term on the right is the KK term, the second term is 
the high-energy modification. The other modification is a straightforward high-energy correction of the background 
quantities H and p via the modified Friedmann equations. 

In the 1+3-covariant approach, perturbative quantities are projected vectors = 0) and projected symmetric 

tracefree tensors, 



(ob> 



ha'^hb'^ — -habh'^'^ 



W,d- (34) 



These are decomposed into (3D) scalar, vector and tensor modes as [6] 

Va = IiaV + Va, (35) 

Wab = D(„Db)W^ + D^aWb) + Wab , (36) 

where an overbar denotes a (3D) transverse quantity (D"t4 = = T>^Wab)- Purely scalar perturbations are charac- 
terized by 

Va = Wa = Wab = 0, (37) 

and scalar quantities are formed via the (3D) Laplacian: V = D'^DaV = D^V. Purely vector perturbations are 
characterized by 

Va = K, Wab = ^{aWb) , curlD„/ = -2/a;„ , (38) 
where Wa is the vorticity, and purely tensor by 

Da/ = = K , Wab = Wab ■ (39) 

The KK energy density produces a scalar mode D^p* (which is present even if = in the background). The KK 
momentum density carries scalar and vector modes, and the KK anisotropic stress carries scalar, vector and tensor 
modes: 

ql = Baq* + q*a , (40) 

<b = D(aD6)7r* + ^^aK) + Kb ■ (41) 



Density perturbations 

We define matter density and expansion (velocity) perturbation scalars, as in general relativity, 



P 



and then define dimensionless KK scalars [11], 



u = -dV* , Q = -B^q* , n = -B\* . 



(42) 



(43) 



A non-adiabatic isocurvature mode is associated with the KK fluctuations. We can see this from the l+S-covariant 
expression which determines the non-adiabatic perturbations in the matter plus KK "fluid" [19] : 



' 3 



- + w; + c: I ^ 



{3pU - 4p*A) , 



(44) 



where Cg = p/p, and the matter perturbations are assumed adiabatic. This mode is in general present, in particular 

when p* = in the background. However the mode will in general decay and be suppressed at low energies. 

This can be related to the curvature perturbation TZ on uniform density surfaces, which is defined in the metric-based 
perturbation theory. The associated gauge-invariant quantity 
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(45) 



may be defined for matter on the brane (^™), for the KK "fluid" (^*) if p* ^ 0, and for the total, efi'ective fluid (C*^^). 
If p* 7^ in the background, 



feff 



r + 



4p* 



(r - r ) • 



(46) 



.3(p + p)(l + p/A)+4p* 

In the case where p* = in the background, the evolution of the total curvature perturbation on large scales is [20]: 



eeff 



P + P 
p + \ 



Sp* 



{p + p){l + p/X) 



(47) 



For adiabatic matter perturbations, ^™ = 0, independent of brane- world modifications to the field equations, since 
this result depends on energy conservation only [27]. However, ^''^^ even for adiabatic matter perturbations. The 
KK effects on the brane contribute a non-adiabatic mode, although ^ '^'^ ^ at low energies. 
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FIG. 1. The evolution of $ along a fundamental world-line for a mode that is well beyond the Hubble horizon at 
N = ln(a/oo) = 0, about 50 e- folds before inflation ends, and remains super-Hubble through the radiation era. A smooth 
transition from inflation to radiation is modelled by w = | [(2 — a) tanh(A'^ — 50) — (1 — a)], where a is a small positive parameter 
(chosen as a = 0.1 in the plot). Labels on the curves indicate the value of po/A, so that the general relativistic solution is the 
dashed curve (po /\ — 0). For po/A 3> 1, inflation ends at A'" = 50 — 2 ln[(l — 2a)/3] ~ 47.4, and at A'" = 50 in general relativity. 
Only the lowest curve still has p/A S> 1 at the start of radiation-domination (A' greater than about 53), and one can see that 
$ is still growing. 



The covariant density perturbation equations on the brane reduce to [19] 



A = 3wHA- {l+w)Z , 

Z = -2HZ - ( D^A - K^pU - Ik'p 



l + {A + '6w)i{ 
A 



( ^ 


-1 




P. 



A, 



(48) 

(49) 



7 



U={3w- l)HU 



= (3 - lw)HQ 



Act 



3a 



U 



HA 
1 



4p* 
4c2 



1 + «; 



Z - aD^Q , 



--3(l+«;)- 



A. 



(50) 
(51) 



It follows that the system closes on super-Hubble scales, since the KK anisotropic stress term 11 occurs only via its 
Laplacian [11]. KK effects introduce two new isocurvature modes on large scales (associated with U and Q), as well 
as modifying the evolution of the adiabatic modes [19,23]. A simple illustration of this modification is shown in Fig. 1 
(from [19]). The variable 



(52) 



is a covariant analogue of the Bardeen metric perturbation variable The large-scale behaviour of $ through the 
inflationary and radiation eras is compared for different values of the inflationary energy scale relative to the brane 
tension. In general relativity, <& is constant in both eras, whereas $ grows during the high-energy regime on the brane. 

Although the density perturbations can be found on super-Hubblc scales, the Sachs- Wolfe effect requires 11 in order 
to translate from density to metric perturbations. In the longitudinal gauge of the metric perturbation formalism, 
the gauge-invariant metric perturbations at last scattering are related by 



— $ 



H 



(53) 



where the radiation anisotropic stress on large scales is neglected, as in general relativity, and Sn* is equivalent to the 
covariant quantity 11. In general relativity, the right hand side is zero. The brane-world corrections to the general 
relativistic direct Sachs- Wolfe effect arc given by [20] 



ST 



Prad 
3 \ Pcdm 



S* - K^a'^Sn* + 



2k' 



j da aJ/'^ 5iT* , 



(54) 



where S* is the KK entropy perturbation (determined by 6p*). The KK term 6ix* cannot be determined by the 4D 
brane equations, so that 5T/T cannot be evaluated on large scales without solving the 5D equations. 

A simple approximation to S-k* on large scales is discussed in [22] and the Sachs- Wolfe effect is estimated as 



5T_ 



Sir* 
P 



2/3 



ln(teq/i4) 



(55) 



where is the 4D Planck time and tin is the time when the KK anisotropic stress is induced on the brane, which is 
expected to be of the order of the 5D Planck time. 



Vector perturbations 



The vorticity propagation equation on the brane is the same as in general relativity, 

Wa + 2ifWa = -^CUrlAa . (56) 

Taking the curl of the conservation equation (24) (for the case of a perfect fluid, 5o = = i^ab), and using the identity 
in Eq. (38), one obtains 

cm\A^ = -QHclwa, (57) 
as in general relativity, so that Eq. (56) becomes 

6ja + [2-icl]Huj^ = Q, (58) 

which expresses the conservation of angular momentum. In general relativity, vector perturbations vanish when the 

vorticity is zero. By contrast, in brane-world cosmology, bulk effects can source vector perturbations even in the 
absence of vorticity [6] . This can be seen via the divergence equation for the magnetic part Hab of the 4D Weyl tensor 
on the brane: 
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4 2 * 1 2 , -* 



(59) 



where Hab = ^(aHb) ■ Even when uja = 0, there is a source for gravimagnetic terms on the brane from the KK quantity 
curlg*. 

We define covariant dimensionless vector perturbation quantities for the vorticity and the KK gravimagnetic term: 



aa = auja, Pa = "Curl q* . 



(60) 



On large scales, we can find a closed system for these vector perturbations on the brane [6]: 



S„ + (1 - Zcl) Haa = 0, 
^„ + (1 - 3w)Hpa = Ih 



4 (3c^ - 1) ^ - 9(1 + w) 



2P 

X 



aa 



(61) 
(62) 



Thus we can solve for and /3a on super-Hubble scales, as for density perturbations. Vorticity in the brane matter is 
a source for the KK vector perturbation 0a on large scales. Vorticity decays unless the matter is ultra-relativistic or 

stifFcr {w > i). and this source term typically provides a decaying mode. There is another pure KK mode, independent 
of vorticity, but this mode decays like vorticity. For w = p/p = const, the solutions are 



Pa = Ca 



3iu-l 



3w-l 



do 



3/90 \ao) 



2(3u)-l) 



(63) 
(64) 



where ba = Q = Ca and e^, = (1, 0) for (w i , ^^; = i). 

Inflation will redshift away the vorticity and the KK mode, which is consistent with the analysis in [18] of vector 
perturbations generated during inflation. 

Tensor perturbations 

The covariant description of tensor modes on the brane is via the shear, which satisfies the wave equation [6] 



D Gab — CTab — 5Haab — 
= {Kb + "^HTtlf) . 



2A+iK2|p_3p_(^ + 3p)£| 



Cab 



(65) 



Unlike the density and vector perturbations, there is no closed system on the brane for large scales. The KK anisotropic 
stress 7f*j is an unavoidable source for tensor modes on the brane. These modes and their effect on the CMB are 
discussed in the following section. 



III. CMB ANISOTROPIES IN THE BRANE- WORLD 



The perturbation equations in the previous section form the basis for an analysis of scalar and tensor CMB 
anisotropics in the brane-world. The full system of equations on the brane, including the Boltzmann equation for 
photons, has been given for scalar [23] and tensor [24] perturbations. But the systems are not closed, as discussed 
above, because of the presence of the KK anisotropic stress tt*;,, which acts a source term. For example, in the tight- 
coupling radiation era, the scalar perturbation equations may be decoupled to give an equation for the gravitational 
potential defined by 



£„6 = D/„Dm$ 



{a^by 



(66) 



In general relativity, this equation in $ has no source term, but in the brane-world there is a source term made up of 
TT*;, and its time-derivatives. At low energies (p <C A), and for a fiat background [K = Q), the equation is [23] 



9 



const 



A 



1 */ / 2 3 1 



(67) 



where x = k/{aH), a prime denotes d/dx, and tt^ are the Fourier modes of $ and tt*;,. In general relativity the 
right hand side is zero, so that the equation may be solved for and then for the remaining perturbative variables, 
which gives the basis for initializing CMB numerical integrations. At high energies, earlier in the radiation era, the 
decoupled equation is fourth order [23]: 



729x^^'l" + 3888a;$fc" + (1782 + 54a;^)$fc + 144x$'fc + (90 + x'^)^k = const 



ini 



243 ( I + 



810 /^7r|Y" ^ 18(135 + 2a:'^) (it^' _ 30(162 + (it^ ^ + 30(162 + x^) ^7r| 



X \p J x^ \P J X \P 



(68) 



The formalism and machinery are ready to compute the temperature and polarization anisotropics in brane-world 
cosmology, once a solution, or at least an approximation, is given for 7r*f,. The resulting power spectra will reveal the 
nature of the branc-world imprint on CMB anisotropics, and would in principle provide a moans of constraining or 
possibly falsifying the brane-world models. Once this is achieved, the implications for the fundamental underlying 
theory, i.e. M theory, would need to be explored. 

However, the first step required is the solution or estimate of tt*;,. This solution will be of the form, expressed in 
Fourier modes (and assuming no incoming 5D gravitational waves): 

7r^(i)cx J dt' g{t,t')Fk{t') , (69) 

where Q \s a retarded Green's function evaluated on the brane. The functional will be determined by the covariant 
brane perturbation quantities and their derivatives. It is known in the case of a Minkowski backgroimd [28], but not 
in the cosmological case. Once Q and Fi~ are determined or estimated, the numerical integration in Eq. (69) can in 
principle be incorporated into a modified version of a CMB numerical code. 

In order to make some progress towards understanding brane-world signatures on CMB anisotropics, we need to 
consider approximations to the solution. The nonlocal nature of tt^, as reflected in Eq. (69), is fundamental, but 
is also the source of the great complexity of the problem. The lowest level approximation to tt^ is local. Despite 
removing the key aspect of the KK anisotropic stress, we can get a feci for its influence on the CMB if we capture at 
least part of its qualitative properties. The key qualitative feature is that inhomogeneity and anisotropy on the brane 
are a source for KK modes in the bulk which "backreact" or "feed back" onto the brane. 

The simplest case is that of tensor perturbations. The transverse traceless part of inhomogeneity and anisotropy 
on the brane is given by the transverse traceless anisotropic stresses in the geometry, i.e. by the photon anisotropic 
stress TTab and the shear anisotropy aab- The photon anisotropic stress in turn is sourced by the shear to lowest 
order (neglecting the role of the octupole and higher Boltzmann moments), so that Fk ~ F[ak\. The simplest local 
approximation which reflects the essential qualitative feature of the spin-2 KK modes is [24] 

K^TfJ, = -CHaab , (70) 

where C is a dimensionless KK parameter, which is assumed to be a comoving constant in a first approximation, 

C = 0. (71) 

The limit C = corresponds to no KK effects on the brane, and C — = gives the general relativity limit. 

For tensor perturbations, there is no freedom over the choice of frame (i.e. u"), and thus there is no gauge ambiguity 
in Eq. (70). However, for scalar (or vector) perturbations, this relation could only hold in one frame, since tt*;, is 
frame-invariant in linear theory while Gab is not: 

U^^U^ + V" <b ^ , aab ^ (Tab + ^{aVb) ■ (72) 

Thus for scalar perturbations, we would need an alternative, frame-invariant, local approximation. 

The approximation in Eq. (70) has the qualitative form of a shear viscosity, which suggests that KK effects lead to 
a damping of tensor anisotropics. This is indeed consistent with the conversion of part of the zero-mode at Hubble 
re-entry into massive KK modes [16,17]. The conversion may be understood equivalently as the emission of KK 
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gravitons into the bulk. This leads to a loss of energy in the 4D graviton modes on the brane, i.e. to an effective 
damping. The approximation in Eq. (70), although local, therefore also incorporates this key feature qualitatively. 

The 1+3-covariant transverse traceless quantities are the electric {Eat) and magnetic parts of the brane Weyl 
tensor, the shear, and the anisotropic stresses. They are expanded in electric (Qab) and magnetic (Qit^) parity tensor 
harmonics [29], with dimensionless coefficients: 



E, 



ab 



= E 



Hab 



f^ab 



(k) 



k 



HkQ^^^ + HkQ':,' 



(fe) 



k 

TTab = pE [^kQ'^ab + ^kQ^ab 
k 

<b = pY.ViQ^^+nQ^^ 



(73) 

(74) 
(75) 
(76) 
(77) 



Using Hab = curl (Tab, we arrive at the coupled equations [24] 



a2 



{a'k+n(Jk) + ^Ek 



-pT^k 



{Ei + nEk)-k 



3K 



-p^k- 



= "tIy ["^^^l + '^^P^'"' - ^^P' + ^P')P''k - 3(3w + l)p {pT^'k + p'T^k) - 9(3w + l)p'W7rfe] 
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2 

-kK^P*ak - — [pTT*k + (p' + Up) TTfc] , 



(78) 



(79) 



where r is conformal time, a prime denotes d/dr, TL = a' /a, and the equation-of-state parameter w is not assumed con- 
stant. Equations (78) and (79), with all brane- world terms on the right-hand sides, determine the tensor anisotropics 
in the CMB, once TTfe and tt^ are given. The former is determined by the Boltzmann equation in the usual way [29], 
since high-energy corrections are negligible at and after nucleosynthesis. The latter is given by the approximation 
Eq. (70), which gives 



2 * /-nj ^ 

1^ PT^k = -CW^O-fe . 



(80) 



We will also assume K = = p* in the background. The KK parameter ( (together with the brane tension A) then 

controls brane-world effects on the tensor CMB anisotropics in this simple local approximation. 

Note that the 4D metric perturbation variable, Ht, which characterizes the amplitude of 4D gravitational waves, 
is related in flat models to the covariant variables by 



Hxk 



k 



+ 2Ek. 



If the photon anisotropic stress TTfc can be neglected, Eq. (78) implies 

kHrk = -<7fc - (C + 2)nak • 



(81) 



(82) 



IV. CMB TENSOR POWER SPECTRA 
In the tight coupling regime we can neglect the photon anisotropic stress (i.e. Wk =0), and the variable 



Uk = a^+^/^cTfe 



(83) 
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satisfies the equation of motion 



a-i-f/2 

by Eq. (80). In flat models {K = 0) on large scales, the solution is 



uk=0, (84) 



a, = ^fca-(2+C) + B,a-(2+C) J dr a(r)2+^ , (85) 

where A), (decaying mode) and Bk are constants of integration. If we let C — * 0, we recover the general relativity 
solution. 

We can solve Eq. (84) on all scales in the high-energy (p A and a cx r^/^) and low-energy (p <^ X and a oc t) 
radiation-dominated regimes, and during matter- domination (a a r^). The solutions are 



Ukir) = Vkr [ci J(5+^)/6(fcT) -|- C2Y(^5_^_Q/e{kT)] (high energy radiation), (86) 
Ufe(r) = Vkr [c3 J(3+f)/2(fcT) -|- 041^(3+0/2 (^t)] (low energy radiation), (87) 
Ufe(r) = Vkr [c5 J(5+2^)/2(fcT) + C6Y(5+2C)/2(fcT)] (matter domination), (88) 

where Cj are integration constants and Jn^Yn are Bessel functions. The solutions for the electric part of the brane 
Wcyl tensor can be found from Eq. (78) . For modes of cosmological interest the wavelength is well outside the Hubble 
radius at the transition from the high energy regime to the low energy. It follows that the regular solution (labelled 
by ci) in the high-energy regime will only excite the regular solution (03) in the low-energy, radiation-dominated era. 
Performing a series expansion, we arrive at the appropriate initial conditions for large-scale modes in the low-energy 
radiation era: 

+ 2(3 + 0(5 + ^''^ 
_ (4 + (fcr)2(8 + ,^..,.41 

- 2(3TO 4(3 + 0(5 + + ^ ^' ^ ^ 

In the limit C ^ 0, we recover the general relativity results [29]. 

For modes that are super-Hubble at matter-radiation equality (i.e. kr^.^ <^ 1), the above solution joins smoothly 
onto the regular solution labelled by C5 in Eq. (88). For fcreq » 1, the shear during matter domination takes the form 

= _2(3+2f)/2p[(5 + 20/2](fcr)-(3+2f)/V(5+2C)/2(fcT). (92) 

In the opposite limit, the wavelength is well inside the Hubble radius at matter-radiation equality. The asymptotic 
form of the shear in matter domination is then 



r[(3 + 0/2] ^2roqV+'^/' 
fft ~ -j= 



(fc^)-(i+C/2)sin(^fcr- Jc) ■ (93) 
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FIG. 2. The temperature power spectrum for tensor perturbations in brane-world models, using the approximation in 
Eq. (70), with the dimensionless KK parameter. Models are shown with ( = 0.0, 0.1, 0.25, 0.5 and 1.0. The initial tensor 
power spectrum is scale invariant. The background cosmology is the (concordance) spatially flat ACDM model with density 
parameters = 0.035, flc = 0.315, Q,a ~ 0.65, no massive neutrinos, and the Hubble constant Hq — 65 kms~^Mpc~^. 

The initial conditions Eqs. (89)-(91), are used in a modified version of the CAMB code [30] to obtain the tensor 
temperature and polarization power spectra [24]. The temperature and electric and magnetic polarization spectra 
are shown in Figs. 2-4 for a scale-invariant initial power spectrum. The normalization is set by the initial power in 
the gravity wave background. Within the local approximation to tt^, the power spectra are insensitive to high-energy 
effects: the ( — curve in Fig. 2 is indistinguishable from that of the general relativity model. For the computations, 
the lowest value of the brane tension A, consistent with the limit A > (100 GeV)*, is used, but the results are largely 
insensitive to the value of A within the local approximation. 

There are three notable effects arising from our approximation to the KK stress: 

(1) the power on large scales reduces with increasing KK parameter C; 

(2) features in the spectrum shift to smaller angular scales with increasing 

(3) the power falls off more rapidly on small scales as C increases. 

Neglecting scattering effects, the shear is the only source of linear tensor anisotropies [29]. For 1 <C ? < 60 the 
dominant modes to contribute to the temperature C;s are those whose wavelengths subtend an angle ~ l/l when 
the shear first peaks (around the time of Hubble crossing). The small suppression in the C;s on large scales with 
increasing ^ arises from the reduction in the peak amplitude of the shear at Hubble entry [see Eq. (92)], qualitatively 
interpreted as the loss of energy in the 4D graviton modes to 5D KK modes. 

Increasing ( also has the effect of adding a small positive phase shift to the oscillations in the shear on sub-Hubble 
scales, as shown e.g. by Eq. (93). The delay in the time at which the shear first peaks leads to a small increase in the 
maximum / for which l{l + 1)C; is approximately constant, as is apparent in Fig. 2. The phase shift of the subsequent 
peaks in the shear has the effect of shifting the peaks in the tensor C;s to the right. For I > 60 the main contribution 
to the tensor anisotropies at a given scale is localized near last scattering and comes from modes with wavenumber 
k I/tq, where tq is the present conformal time. On these scales the gravity waves have already entered the Hubble 
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radius at last scattering. Such modes are undergoing adiabatic damping by the expansion and this results in the 
sharp decrease in the anisotropies on small scales. Increasing the KK parameter ^ effectively produces more damping 
and hence a sharper fall off of power. The transition to a slower fall off in the C;s at I ~ 200 is due to the weaker 
dependence of the shear amplitude on wave-number at last scattering for modes that have entered the Hubble radius 
during radiation domination. The asymptotic expansion of Eq. (92) gives the shear amplitude oc at fixed r, 

whereas for modes that were sub-Hubble at matter-radiation equality Eq. (93) gives the amplitude oc 




Similar comments apply to the tensor polarization and CP , shown in Figs. 3 and 4. As with the temperature 
anisotropies, there is the same shifting of features to the right and increase in damping on small scales. Since 
polarization is only generated at last scattering (except for the feature at very low I that arises from scattering at 
reionization, optical depth 0.03) the large-scale polarization is suppressed, since the shear (and hence the temperature 
quadrupole at last scattering) is small for super-Hubble modes. In matter domination the large-scale shear is (Tfc = 
— fcr/(5 -|- 2C); the reduction in the magnitude of the shear with increasing KK parameter C is clearly visible in the 
large-angle polarization. 
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FIG. 4. The magnetic polarization power spectrum for tensor perturbations for the same brane-world models as in Fig. 2. 



V. CONCLUSION 

In principle, observations can constrain the KK parameter C, which determines the brane-world effect on tensor 
anisotropics in the CMB by controlling the generation of 5D modes within a simplified local approximation, Eq. (70). 
The other brane-world parameter A, the brane tension, is not constrained within this approximation. This may 
indicate that the approximation entails a hidden low-energy assumption, or it may simply be accidental. 

The local approximation to 7r|j introduced in [24] is a first step towards the calculation of the brane-world imprint 
on CMB anisotropics on small scales, starting with the simplest case of tensor anisotropies. In practice, the tensor 
power spectra are not measured, and the prospect of useful data is still far off. What is more important is the 
theoretical task of improving on the simplified local approximation in Eq. (70). This approximation encodes aspects 
of the qualitative features of brane-world tensor anisotropies, primarily the loss of energy in 4D graviton modes via 5D 
graviton emission, which may be expected to survive in modified form within more realistic approximations. However, 
a proper understanding of brane-world effects must incorporate the nonlocal nature of the KK graviton modes, as 
refiected in the general form of Eq. (69). 

Furthermore, it is the scalar anisotropies which dominate the measured power spectra, and it is therefore of even 
greater importance to develop the scalar analysis. In this case, even the first step of a local approximation has 
to confront the problem of the frame ambiguity in Eq. (70). An alternative frame- invariant local approximation is 
needed, as a first step. Nonlocal approximations, and ultimately numerical integration of the nonlocal equation for 
TT^, must be developed. 

Once these more realistic nonlocal approximations are developed and the brane-world imprint on the CMB is 
computed, there are two key tasks: 

• generalize the CMB brane-world computations to include a bulk scalar field and other extensions of the RSII-type 
models; 
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• determine the implications of the branc-world imprint for the fundamental theory, M theory, which underlies 
the key aspects of phenomenological brane- world cosmologies. 
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